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A method is presented for measuring gaseous diffusion coefficients using a gas chromatographic 
technique. Diffusion coefficients were measured for the systems: argon, krypton, oxygen, and nitrogen 
diffusing into helium at temperatures from 77 to 400 K. 
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1 . Introduction 

Accurate knowledge of molecular diffusion coeffi- 
cients and their temperature dependence is important 
in obtaining fundamental information on forces be- 
tween unlike molecules. Until recently most diffusion 
measurements were conducted over short temperature 
intervals because the conventional methods such as 
the Loschmidt or the Stefan, require the thermostating 
of large, static systems. Walker and Westenberg [l] 1 
have developed a method which involves measuring 
concentration profiles downstream of a point source of 
tracer gas in a uniform laminar jet of carrier gas. 
Diffusion coefficient for permanent gas pairs over a 
temperature interval 300 to 1100 K have been measured 
using this technique. 

Determining diffusion coefficients using gas cro- 
matographic techniques is relatively new. Taylor [2] 
derived an expression for the radial diffusion term for 
a liquid dispersion flowing in a tube. The same expres- 
sion was derived by Giddings [3] and Golay [4] as a 
special case of a general theory of gas chromatography. 
Giddings and Seager [5] have measured diffusion 
coefficients using this method. 

The method reported in this paper is also based on 
gas chromatography techniques. In brief, a trace 
amount of one gas was injected into a carrier gas 
stream of different composition flowing through a 
column at a constant flow rate. The concentration 
profile of the gas dispersion was determined at two 
points downstream from the point of sample injection 
by thermal conductivity detectors. The binary diffu- 
sion coefficients were calculated from the dimensions 
of the recorded peaks and time required for the peak 



maximum to travel the distance between the two 
detectors. 

Diffusion coefficients were measured for the sys- 
tems: argon, krypton, oxygen and nitrogen diffusing 
into helium at temperatures from 77 to 400 K. 



2. Theoretical 



2.1. Fundamental Relations 

As has been shown by theoretical studies [2, 3, 4] 
when a nonuniform gaseous dispersion flows through 
a cylindrical tube, the effect of radial diffusion acting 
in conjunction with the parabolic velocity profile is 
accounted for by an effective longitudinal diffusion 
coefficient a 2 , given by 
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where v is the mean flow velocity, a is the inside radius 
of the tube, and D is the diffusion coefficient. If c 
represents the mean concentration in a cross section 
at distance x from the moving centroid of the disper- 
sion (x is taken positive downstream from the centroid) 
then c satisfies the one-dimensional diffusion equation 
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1 Figure in brackets indicate the literature references at the end of this paper. 



First let us assume an ideal system in which the 
experiment is initiated by instantaneous injection at 
time t = of n moles of diffusing substance into the 
carrier gas stream at a fixed injection point, and 
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further assume that the external parameters affecting 
the flow (e.g., temperature, pressure, and tube radius) 
are all constant. Equation (2) applies after a time which 
is long compared with a 2 /D. Its solution appropriate 
to these initial conditions is given by 
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Equation (3) defines a concentration distribution 
which is instantaneously Gaussian and spreading 
with time. In the actual system the assumption of 
constant flow parameters applies only to the diffusion 
column proper. However, the experimental system 
was designed such that the concentration distribution 
was always instantaneously Gaussian within the dif- 
fusion column. Equation (3) can be applied under 
these conditions if the zero of time is redefined. 
Instead of taking £ = as the instant of sample injec- 
tion, it must be so defined that eq (3) is consistent 
with an analysis of the peaks recorded from two or 
more separate detectors along the tube. 

Equation (3) expresses concentration in terms of a 
frame of reference that moves with the mean velocity 
of the carrier gas stream. In order to express concentra- 
tion as a function of time at a fixed detector, we must 
transform eq (3) by means of the relation 



x = v(t -t). 



(4) 



where to is that time at which the centroid of the 
dispersion arrives at the detector. When we make 
this substitution, the concentration at a fixed detector 
is given by 
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2.2. Peak Analysis 

Assuming linear detector response, the height h of a 
recorder trace of the signal at time t is proportional to c 
as expressed by eq (5). If ho is the height at time to, it 
follows from eq (5) that 
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From this equation it is seen that a straight line with 
the slope v 2 /4>a 2 should be obtained when 
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is plotted versus (to — t) 2 . Such a plot is useful for 
demonstrating that the peak shape conforms to theory. 
However, for determining the diffusion coefficient an 
approximation based on eq (6) is easier to apply while 
providing sufficient accuracy. 
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Since e is small (typically, e 2 < 10~ :i ), eq (8) can be 
expanded in powers of e. Neglecting powers higher 
than the second, 

(l + e'W-e^ln^'-l^-ln^O. 

This equation is quadratic in r\. Let A77 be the differ- 
ence of its two roots. Expanding In -7- as a power series 

in e 2 , the solution of eq (8), correct through terms in 
e 2 , can be written 

4 1ny= (Ar ? ) 2 + e 2 [2(A7)) 2 -i(AT ? ) 4 -]]. (9) 

Taking typical values for e and A77 (typically, A^ 2 < 10), 
it can be readily seen that all but the first term on the 
right side of this equation can be neglected when 
accuracy to a few parts per thousand is adequate. 
If we let At represent the width of the recorder trace 
at height h, we can make use of eqs (7a) and (7b) to 
write 



At/ = 



A; 



aVTo 2 



Substituting in eq (9) and retaining only the first term 
on the right we obtain the approximation 



to In 



h ~ 4a 2 V 2 



If £0 In -7- is plotted versus I — 



it is seen that a 



straight line with the slope v 2 /ia 2 should be obtained. 
This plot should superimpose accurately on that based 
on eq (6). It is obvious that eq (10) is easier to apply in 
determining a. 

Since to, the time of arrival of the peak at a detector, 
appears on the left in eq (10), it is still necessary to 
define the zero of time. The difference of the to values 
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This approximation can be derived by introducing 
the characteristic diffusion distance s defined by 

s = aVu } , (7a) 

the dimensionless variable 17 defined by 

v(t-to) = 2r]s, (7b) 

and the small dimensionless parameter e defined by 

s = €vt . (7c) 

When eq (6) is rewritten in terms of e and r/, the result is 



(10) 



at two detectors along the tube can be measured di- 
rectly. Their ratio can be obtained by applying eq (10) 
to each of the detectors. Solving for the two to from 
their difference and their ratio is equivalent to deter- 
mining the effective zero of time. 

Another method which, although tedious, yields a 
with good precision because random fluctuations in 
the recorder signal are smoothed by integration, is 
based on the ratio of peak area to peak height. The 
equation 
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f equivalent to eq (6), can be reduced to a standard 
integral form (7), by the substitution y 2 = v 2 t/4<(x 2 , with 
the result 
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2.3. Experimental Procedure 



In figure 1 is shown a schematic drawing of the 
apparatus. A 150-ft emptyV4-in O.D. copper column 
was divided into three equal segments; a thermostatic 
column, Ci, the diffusion column, C 2 , and a third 
column, C;j, which allowed the carrier gas to come to 
room themperature before its flow rate was measured 
by a soap-bubble flowmeter, F. The thermostating 
and diffusion columns were immersed in liquid nitro- 
gen for the low temperature measurements and an oil 
thermostat, controlled to ±0.02 °C, for the other 
measurements. At the beginning and end of the diffu- 
sion column were thermal conductivity detectors, 
Di and D 2 , located in line with the incoming gas flow. 
The detectors, shown in figure 2, were constructed 
from V4-in Tee connectors for tubing in which one arm 
was cut off to allow the sensing element, a Gow Mac 
type 9225 filament support (see footnote 2) to be sol- 
dered into position so that the sensing element was on 
the central axis of the diffusion column. The two 



r^Oir^ — ft 





n-AWWW\rt^ 



<=t= 



Figure 2. Detector. 



detectors formed adjacent arms of a Wheatstone bridge 
circuit. The resultant voltages were plotted as a func- 
tion of time using a 1 mV full scale recorder. 

The flow of helium carrier gas from a high pressure 
tank, H, was controlled by a Moore [6] 2 constant- 
differential type flow controller, M, a needle valve, 
Ni, upstream from the diffusion column and two 
needle valves, N 2 and N3, downstream from the dif- 
fusion column. An 0.1 ml of sample gas from a high 
pressure tank, S, was injected into the thermostating 
column by means of a F & M [8] gas injection valve, 
J. Commercially available gases with a stated purity 
of better than 99 mole percent were dried by passing 
them through a molecular sieve column (Ti and T 2 ). 
The column pressure was read from a Wallace and 
Tiernan [9] absolute aneroid barometer, G. 

The effective radius, a, was calculated at room 
temperature by equating the mass flow rate measured 
by the soap-bubble meter to the mass flow rate in the 
diffusion column, 7rP c a 2 L/t, where P c is the pressure 
of the carrier gas in atmospheres and t is the time 
required for an air peak to travel the distance, L, be- 
tween the two detectors. The value of a agreed within 
one part in two hundred for all the measurements 
taken. 

For the high and low temperature measurements 
the mass flow rate in the diffusion column was used to 
determine the temperature of the carrier gas, the 
column's dimensions being corrected for thermal 
expansion. If the carrier gas was not at the tempera- 
ture of the thermostat the flow rate was reduced until 
it was. 

A preliminary diffusion measurement was made for 
each system to determine peak heights. The attenua- 
tion was then adjusted by means of a 10-turn, 10,000 ft 
Heli-pot so that for all the following measurements 
the peak heights would be the maximum height al- 
lowed by the recorder chart paper. Duplicate runs were 
made by injecting successive amounts of sample gas 
into the carrier gas stream, the injections being at a 
time interval such that no overlapping of peaks oc- 
curred. At liquid nitrogen temperature, it was possible 
for as many as four peaks to be in the diffusion column 
simultaneously. This enables four velocity measure- 
ments to be made during the course of a run. The time 



FIGURE 1. Schematic diagram of flow system for measuring diffusion 
coefficients. 



2 Certain commercial materials and equipment are identified in this paper in order to 
specify adequately the experimental procedure. In no case does such identification imply 
recommendation or endorsement by the National Bureau of Standards, nor does it imply 
that the material or equipment identified is necessarily the best available for the purpose. 
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between peaks was checked with a stop watch. The 
heights and widths of the peaks were measured to a 
precision of one part in two hundred fifty with an 
engineer's scale. 

The quantity a 2 was calculated from the expression 



a 



il- 1 
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(13) 



where S 2 and Si are the slopes of the In Hoi h — (At) 2 
plots and (£02 — £01) is the time required for the peak 
maximum to travel the distance between detectors. 
The quantity a 2 was also calculated from 
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where A 2 and A x are the areas under the second and 
first peak respectively and v the velocity of the carrier 
gas in the diffusion column. Equation (14) which follows 
from eq (11), gave the same values for a 2 as eq (13). 

3. Results and Discussion 

In figure 3 are shown plots of eqs (6) and (10) for one 
experiment on the diffusion of nitrogen in helium at 
30 psia and 25°. The filled squares represent values 
calculated from eq (10) for the first peak with £ i equal 
to 23.59 min and the open squares are values obtained 
from the second peak with £02 equal to 48.81 min. The 
open circles represent values calculated from eq (6) 
from the leading half of the first peak while the filled 
circles are those of the leading half of the second peak. 
The half-filled circles right represent values calcu- 
lated from eq (6) from the trailing half of the first peak 
while the half-filled circles left are the corresponding 
values for the second peak. All experimental points 
fit the theoretical line within a precision of 0.5 percent. 
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There are no indications of skewing or tailing of either 
peak. 

Values for the diffusion coefficients at different 
temperatures are given in table 1. The pressure of the 
carrier gas in the diffusion column varied from 1.1 to 
1.6 atm for these determinations. These values have 
been corrected to standard atmospheric pressure. 
Each diffusion coefficient was obtained as an average 
of eight separate determinations. The values of 0.732 
and 0.688 cm 2 /s for the diffusion of argon in helium and 
nitrogen in helium at 25 °C respectively compare 
favorably with values of 0.725 and 0.688 cm 2 /s obtained 
by Walker and Westenberg [1] for the same gas pairs. 

Table 1. Effect of Temperature on the Diffusion Coefficient 



Gas pair 


Temp. K 


Z)(cm 2 /s) 


He-A 


77.2 


0.0710 ±0.0009* 


He-A 


298 


.7335± 


.004 


He-A 


334 


.8890± 


.006 


He-A 


357 


.9917 ± 


.006 


He-N 2 


77.2 


.0725 ± 


.001 


He-N 2 


296 


.678 + 


.005 


He-N 2 


321 


.8171 ± 


.005 


He-N 2 


348 


.9251 ± 


.007 


He-N 2 


370 


1.030 ± 


.006 


He-Kr 


298 


.6491 ± 


.004 


He-Kr 


322 


.7372 ± 


.004 


He-Kr 


341 


.813 ± 


.005 


He-Kr 


366 


.904 ± 


.006 


He-0 2 


298 


.7361 ± 


.005 


He-O, 


320 


.8472 ± 


.006 


He-0 2 


365 


1.041 ± 


.008 



*90% confidence limits for the mean. 



A least square calculation for the expression 
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FIGURE 3. Plot of eqs (6) and (10) for the diffusion of nitrogen in 
helium at 30 psia and 25 °C. 



where Z) 2 98 is the diffusion coefficient at 298 K, D T is 
the diffusion coefficient at temperature, T, and mis a 
constant. The following values for m were obtained 
with indicated precision being for a 90 percent confi- 
dence limit: He-N 2 , 1.69 ±0.02; He-0 2 , 1.67 ±0.06; 
He-Kr, 1.71 ±0.06 and He-Ar, 1.72 ±0.006. These 
values may be compared with values obtained by 
Seager, Geerston, and Giddings [5]: He-N 2 , 1.75 ±0.12; 
He-0 2 , 1.70 ±0.08 and He-Ar, 1.73 ±0.07. The two sets 
agree within the given experimental errors. In figure 
4 is shown the log diffusion coefficient-log temperature 
plot for argon diffusing into helium. 

There are several advantages to this method. 

(a) The use of two detectors instead of one eliminates 
the need for a zero time correction, since it is not 
necessary to know the shape of the peak at the point 
of sample injection. The only requirements are that 
the concentration profile be Gaussian when the peak 
reaches the first detector, that the thermal conductivity 
detectors do not distrub the concentration profile, 
and that the carrier gas be at thermostat temperature 
since the sample injection value is at room temperature. 

(b) Diffusion coefficients over a wide range of tem- 
perature may be measured since the apparatus may be 
used in a liquid nitrogen bath or in an oven without 
any change in design. 
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FIGURE 4. Plot of log diffusion coefficient-log temperature for argon 
diffusing into helium. 

(c) Time measurements may be made accurately. 
Since the velocity of the carrier gas appears as a square 



term in the diffusion equation it is important to be 
able to measure this quantity accurately. This can be 
done by measuring the time required for the peak 
maxima to travel the measured distance between the 
two detectors. By injecting more than one peak into 
the diffusion column it is possible to determine whe- 
ther the velocity of the carrier gas has changed during 
a measurement. 

(d) Deviations of experiment from the assumed 
model may be determined. Equation (6) provides a 
means for analyzing the peak for any disturbance due 
to gas absorption on the walls of the diffusing column 
or for any disturbance to the concentration profile 
due to the gas injection and detection systems. 

It is conluded that this is a simple, rapid, widely 
applicable method for the determination of gaseous 
diffusion coefficients. 
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